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Abstract 

A many-body wave function is approximated by a product of two func- 
tions: the wave function (f> depending on the particle coordinates and the 
function x depending only on the value of interparticle interaction poten- 
tial. For the given (j) an ordinary linear differential equation for x is derived 
by averaging the Hamiltonian over the constant interparticle interaction 
potential surface. Generalized Hartree-Fock equations containing correla- 
tion effects are obtained. To test the proposed technique the ground 1^5o 
and excited 2''5'i states of two-electron ions from H~ up to Ne*"*" are cal- 
culated. In all cases the calculated energies are more accurate than those 
obtained with the Hartree-Fock theory even taking as (j) the symmetrized 
product of electron wave functions in the Coulomb field of nucleus com- 
plitly disregarded the electron-electron interaction. Variation of factors 
in the one-particle wave function exponents leads to the results close to 
those of configuration interaction approach. 

PACS: 31.15.-p, 31.15.Ar, 31.25.Eb, 31.25.Jf 

1. INTRODUCTION 

Direct solution of the many-body Schrodinger equation has been performed 
only for two electron ions [1,2]. For more complicated cases dimensionality of 
the equation and impossibility to separate variables in it force to search for ap- 
proximating solutions. There are two main approaches to find these solutions: 
the Hartree-Fock approximation (HF) [3,4], and the density functional theory 
(DFT) [5,6]. The first one is based on the approximation which considers elec- 
tron moving in an average field of other electrons. The wave function in this 
approximation is presented by a symmetrized product of one particle functions, 
so the probability to find an electron at some space point does not depend on 
the positions of the rest electrons. To improve the Hartre-Fock theory one has 
to introduce this dependence known as the correlation in electron movement 
[7]. It can be achieved through configuration interactions, perturbations or 
coupled cluster theory [8]. DFT rests on the density functional theory where 



exchange and correlation arc introduced by a functional of electron charge den- 
sity. Unfortunately, the exact form of the functional is not known and some 
approximations to it are used. In both approaches one needs to solve non-linear 
one-particle equations using an iteration procedure of the self-consistent field. 
In resent years, powerful computers and sophisticated computational algorithms 
have made a remarkable progress in this area. Nevertheless, the final solution 
of the problem has not been achieved yet, and developing new approaches, per- 
mitting to view this problem from another perspective, is being continued [9]. 
Searching for highly precise functions containing thousands variational param- 
eters as well as simple but accurate functions for heliumlike ions has not been 
stop yet also [10-25]. 

The aim of the present paper is to propose a new approach to solve the 
many-body Schrodingcr equation going beyond HF theory. The idea of this 
work has been stated in [26]. We consider the value of the interparticle in- 
teraction potential as a variable and represent a many-body wave function as 
a product of a function (f> depending on particle coordinates and a function x 
depending only on the value of interpaticle interaction potential. Function x 
due to its dependence on this potential introduces correlation in the many-body 
wave function omitted in 4>. For a given (f>, using averaging the many-particle 
Hamiltonian over constant interaction potential hypersurface, the Schrodinger 
equation is transformed into a one-dimensional linear differential equation for %. 
In its turn, for any x one can find the best (f> which minimizes the total energy 
of the many-body system. Below, we develop the correspondent equations and 
apply them to calculate energies of two-electron ions. 



The Schrodinger equation of interacting particles can be written in the form 



where Hq and Vint are contributions to the Hamiltonian from non-interacting 
particles and the interparticle interaction, correspondently. 



r is a vector with components r^, i numerates particles, F is a potential of an 
external field. For the Coulomb interaction 



II. EQUATIONS 
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where nj = \ri - rj\. 

Let us introduce a function of the form 



V'(r) = <j){r)x{p). 



(4) 
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The function </> is supposed to be symmetrized in respect to the particle permuta- 
tions. The type of symmetrization must reflect the total spin of the considering 
system [27]. Substitution of (H]) in (P) leads to 



rl 



dp 



ih{p) + l/p)0(r)x(p) = E<p{v)x{p) (5) 



It is easy to find 



(6) 
(7) 



In ([7]) it was taken into account that the Coulomb potential is satisfied to the 
Laplace equation. 

The demands for function x to minimize the value of the functional 



for a preset 4> leads to the equation for x 

- Up)^-Pf - "b)^ + (Mp) + yp)x{p) = Exip), 



dp^ 



dp 
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S{P) 



and the integration is performed over a constant interaction potential surface. 
The boundary conditions for x follow from the demand for il) to be finite in the 
whole space. 

Eq. ^ is a linear ordinary differential equation for x- The term hip) is the 
contribution of Hq to the total energy when the system is in the state (p. The first 
and the second terms describe additional contributions to kinetic energy which 
doesn't enter hip) and appears due to the interparticle interaction 1/p. If the 
interaction is absent and is an eigenstate of with energy e, the additional 
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contribution must be zero that occur when x = const, h is independent of p 
and equals to e, and E = e. For interacting particles, x can be a constant only 
if 4> would be an exact solution of the Schrodinger equation. The function x 
we will name the correction function. The type of correction in x depends on 
(j). If (j) corresponds to non-interacting particles x contains corrections due to 
the whole interaction, in case of HF (p only correlation part of the interaction, 
in the case of CI the part of correlation not taken into account in the used CI 
approximation. 

Switching on the interaction between particles located on the same constant 
interaction potential surface does not change the particle motion on this surface 
and can accelerate the particles only along the normal to the surface. Eq. 
([9]) takes into account this acceleration. However, the normal to the surface 
at point r can not be parallel to the normal at point r', so the acceleration 
obtained by the particle at r can have a tangential component at r'. Eq. ([9]) 
disregards this effect and relative value in the change of the electron movement in 
normal and tangential directions due to electron-electron interaction for a given 
X will determine the accuracy of the proposed approximation. To improve the 
approximation one can use common techniques choosing the better (j>. 

The best (j) for independent particle approximation can be obtained by min- 
imization of functional ([5]) considering x as a known function. Presenting (j) by 
a determinant of one-particle functions ipi, one can obtain a set of equations 
analogous to the Hartree-Fock ones 

[-^V? + y(ri)](p,(ri) +^T^.iK,(ri)Vi^,(ri) + 

jk 

w^f (riV,(ri)]=£;^,(ri), (11) 

where 

k 

ry(ri) = , u,j{ri) = {(t)i\^x*^kP^\(t>j)i, (12) 

k 

and is the reverse matrix to r. Integration in (|12p is performed over coor- 
dinates of all particles except particle 1, (j>i means that function ipi in (j> is re- 
placed by unity. The set of equation pT|) resembles the Hartree-Fock equations. 
A difference arises from the fact that minimization in respect to one-particle 
functions is performed with weighted function x depending on the coordinates 
of all particles. It leads to the appearance of non-diagonal matrix elements in 
r and replacing the interpaticle Coulomb interaction 1 /ri2 of the Hartree-Fock 
equations by the screened interaction v'^^^. Through this interaction the corre- 
lation in the Hartree-Fock equations is introduced. For x = 1, equations (llip 
are common Hartree-Fock equations. 

Another possibility to improve cj) avoiding the need to solve HF equation 
occurs for (j) depending on some parameters. In this case eigenvalues of ^ can 
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be considered as the functional on these parameters and their values can be find 
minimizing the functional. This technique will be used in section 3. 

Constructing interaction potential surface, it is useful to note that the form 
of the interparticle interaction ([3]) allows to carry any set of particle positions r 
belonging to constant interpaticle interaction potential surface 1/p to another 
set r' belonging to surface \/p' by coordinate scaling with factor p'/p. It makes 
the generating points for any surface an easy task and restricts itself for initial 
generation by the surface, for example, p = 1. To obtain the point for any other 
surface one can use coordinate scaling again. Besides, the construction of the 
constant interaction potential surface can be presented in a recurrent form. For 
this purpose let us introduce notations 



1 ""^ 1 



where n is the number of particles, and express the interaction potential as 
contributions from n — 1 particles and interaction of the nth particle with the 
rest 

+ (14) 

Pn Pn-1 Qn 

Then integration over the constant interaction potential surface S{pn) for n 
particles can be presented in the form 



dRn = / dKn-i / dr„ (15) 

S(p„) JS(p„_i) Js(q„) 

where c?R„ = dvi . . . c?r„. For a given Rn-i = {ri, . . . , fn-i} integration over 
r„ is performed over the surface determined in 3-dimension space by (|14p. For 
two particles it is the sphere of radius p around position of the first particle, for 
three particles it is a circle formed by crossing the sphere of radius ri3 around 
ri with the sphere of radius r23 around r2, for four particles it is only two points 
at the circle, etc. In addition one can restrict itself by irreducible part of the 
surface results from the identity of the particles. Evidently, with the help of 
particle numeration one can always satisfy the condition 

E-^E- if ^<^' (16) 

' tl ■ ' ik 

and it will be enough to consider only the part of the whole surface, satisfying 
to (Ull). 

As an example, we represent the integral over the constant interaction po- 
tential surface for two particles 



/ rfR2 = / rfdri / dipi / sin6'id6'i 
Js(p) Jo Jo Jo 

/•OO /'27T /'IT 

/ S{ri2 - p)rl2dri2 / dip2 / sin92de2, (17) 

JQ Jo Jo 
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where the integration over the coordinates of the 1st particle is performed over 
the whole space, while the integration over the coordinates of the 2nd particle 
for a given ri is reduced to the integration over a sphere of radius ri2 = p 
around the first particle. 

III. HELIUMLIKE IONS 

As an application of the developed technique we calculated the ground 1^5*0 
and excited 2'^ 5*1 states electronic structure of heliumlikc ions from H~ up to 
Ne^"*". For these ions p = ri2 and ^ is reduced to the form used in a number 
of papers [1,2,17,22,23,25]. The difference between those and the present paper 
consists in the different ways of the wave function calculation. 
We represent 



0(ri,r2) = [(^i(ri)^2(r2)±<y92(ri)(^i(r2)]/V2(T±52) (18) 

where 5* = {Lpi\ip2), for \^S^ 

(^,(r) =2zf'exp(-Z,r) (19) 

is Is function of the electron in the field of nucleus with the charge Zi. For 2'^ 5*1 
the function of type (|19p was used for Lpi , whereas 

(/?2(r) = Zl'^{l ~ Z2r/2) exp (-Z2r/2)/\/2 (20) 

is the 2s function of the electron in the field of nucleus with the charge Z2. 
Factor Zi in function (fTO|) or (|20p is known as an efhcient nuclear charge seeing 
by the electron in the corresponding state, Z without subscript we save for the 
real nuclear charge. For function (|18p Eq. ([5]) can be written in the form 



dp'^ 



2pri dri 2pr2 dri p 

Id'^cj) Id^cp 1 d(f> 1 del) 1 1 1 
2 dri 2 ri dri r2 dr2 ri 7'2 Zp 



dp 



X = E4>x (21) 



Here Z is the nuclear charge, the length unit is Bohr/Z, and atomic units are 
used for other values. 

Integration in (|17p over a constant potential surface, when calculating ex- 
pressions (fTOj) for these case, is reduced to 

/ dR2 = inp / ndn / r2dr2 (22) 

Js(p} Jo J\ri-P\ 

and can be performed analytically. The results are cumbersome and we present 
only their common form 

E ci,fc/e-2^^'' + C2,p'e-^^^+^^-^P + E cs^^p-e-^^^^ (23) 

k I m 
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Eq. ^ with coefficients can be solved analytically for asymptotic cases 
p — !■ and p — > oo. The finite asymptotic solutions are 

2 

Jhnx(p)^x(0)(l + ^ + Y|^+Ob')), (24) 



hm x{p) ^ exp [( ^ - JZ2^„^ - S + So)p] , (25) 

p — ^oo V 

where Zmin — min(Zi, Z2), Eq is energy in state 0. 

Eqs. p4l) and (|25p were taken as the bounder conditions in the numerical 
solution of For the numerical calculations the differential operators were 
approximated with finite differences. Adopting 

dXi ^, Xi+i — Xi-i ^25) 



dp 2Ap 



(27) 



where i numerates points of discrete set {pi\, Xi — x{Pi)i = Pi ^ Pi-i-, Eq. 
([9]) can be represented by n x n tridiagonal matrix with n is the number of 
points in {pi\ 

H^^ = + ^ + hi, (28) 



2Ap2 Zp, 

2Ap2 2Ap 



i^M±l = ± 7^ (29) 



for I < i < n, and 



Hll = 7;^ + — + h^+{-^ + ^]/{l + Ap/2Z + Ap'/UZ-'), (30) 



h 


1 


2Ap2 


Zpi 




-i- 



) exp [( - Z„,„ - Jzl^^-E + Eo)Ap] (31) 



2Ap2 2Ap/ 

in accordance with ([M)) . ([^5]) . Because boundary condition (|3ip depends on 
calculating eigenvalue, it is satisfied in the iteration process. Eigenvalues and 
eigenfunctions of ([9|) were approximated by eigenvalues and eigenvectors of ma- 
trix H in interval < p < 20/Z with 200 points at the interval. The value 
of interval and the number of points guarantee all significant numbers in the 
results presented below. 

For I^Sq states the calculations were performed for three sets of function 
(fT9|) : al) Zi = Z2 = Z; h) Zi = Z2, however, the value is chosen from energy 
minimum; c) Zi ^ Z2 and each the value is chosen from energy minimum. For 
2^51 only the first and the third cases were considered. The calculated energies 
of these states together with those of HF theory for I^Sq, and configuration 
interaction for I^Sq and 2^S'i are shown in Table I. 

As it is seen in Table I in all cases the energies of 1^5*0 states calculated with 
Eq. ([9]) are less than the energies calculated with HF theory. Even for functions 
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(j) completely disregard the interaction between electrons the correction function 
X capable to take into account electron-electron interaction more accurate than 
it is done in HF approximation. However, the achieved accuracy is still far 
from the accuracy obtained by configuration interaction calculations [29] or 
calculations which explicitly introduce the distance between electrons in the 
wave function [1,2,23,25]. It is due to correction function x describe the electron 
motion only between different electron-electron potential surface keeping the 
motion on a surface determined by function (p. 

The results can be improved by choosing the best (j) solving Eq. pT|) . How- 
ever, we attempt to use another way. Function depends on two parameters, 
Zi and ^2, so eigenvalues of ^ are the functionals of these parameters. The 
best one-particle functions of the form (fT9|) and ((20)) minimize the functional. 
It permits us to avoid the solution of non-linear HF equations and makes tech- 
nique based on Q self-sufficient. Table I contains obtained results. We observe 
a slight improvement when both electrons are seeing an equal nucleus electric 
charge. Considerable improvement is observed when different electrons are see- 
ing different effective nucleus charges. These charges and their dependence on 
Z are presented in Fig. 1. 

As for 1^5*0 the results for 2^51 obtained with Zi = Z2 — Z can be sig- 
nificantly improved by introducing different effective nuclear charges for the 
electrons. As it is seeing in Table I, in this case the results are close to those, 
obtained with configuration interaction. The charges and their dependence on 
Z are presented in Fig. 2. 

The function (j) obtained by minimization of E{Zi, Z2) can be considered as 
HF function. Thus, the correction function x can be regarded as correlation 
functions. These functions, normalized so that integral x{p)sip)xip) over inter- 
val < p < 20 equals to 1, are presented in Fig. 3. The functions reduce the 
probability of electrons to be at a short distances showing monotonic growth 
with the distance between the electrons. The behavior of x for 3.11 ions is simi- 
lar, however, the functions become smoother with the growth of Z and tend to 
some limit. The deviation x from the unity for 2^ Si states is less then for 1^5*0 
state suggesting a decrease of correlation effects. It is known to occur due to the 
exchange hole appearing for antisymmetric (j). The function x changes this hole, 
with this change being more significant at electron-electron distances p < 1/Z 
reducing further the small probability to find electrons at these distances. 

The function s{p) determined in ([T0|) is the distribution function of ri2 for 
wave function 0, and function x^X is the distribution function for (j)X- These 
function for l"'^5'o states for H~ and 7Ve*+ are shown in Fig. 4. Besides, Fig. 4 
shows the corresponding correlation holes determined as s(x^ ~ !)■ In general, 
the form of the correlation holes is the same as in [30]. However, the depths of the 
holes decreases from H~ to iVe^+. The distribution functions and correlations 
holes for 2'^ Si of He and A^e^+ are shown in Fig. 5. As for 1^5*0 states the 
depths of the holes decreases with nuclear charge growth and the holes become 
more shallow in comparison with 1^5*0. 

Until now we have considered only the lowest eigenstates of However, for 
given the </> we can calculate higher eigenstates for x- We will demonstrate how 
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these states correspond to real excited states on the example of non-interacting 
electron for which the exact results have been available. In Fig. 6 functions x 
for the lowest three solutions of (O for non-interacting electrons are shown. As 
it is seen, higher eigenfunctions exhibit nodes due to which function (j) based 
on Is functions is transformed in function (j) similar to the function based on 
2s or 3s functions. However, the energies corresponding to these function are 
E = —0.440 and —0.245 instead of exact energies -0.625 and -0.555 a.u. The 
source of pure accuracy is the conservation Is function in computations (jlOp . 
The use of proper function in pI7|) makes it possible to obtain the exact results 
for the excited states. As such, we can use higher eigenstates of (O only for 
a crude estimation of exited states energies and use proper (j) to obtain more 
accurate results. 

IV. CONCLUSIONS 

Eq. ([9]) reduces the many-body Schrodinger equation to the 2nd order lin- 
ear ordinary differential equation for the correction function. It describes the 
changes in electron movement between constant interaction potential surfaces 
due to electron-electron interaction. The changes in electron movement along 
the surface are disregarded in 

Solutions of the Q for two electron ions give the total energies more accurate 
than HF theory does. It is true even in the simplest case where wave functions 
of non-interacting electrons are used for function <p. For </> containing only 
two variations parameters, efficient nuclear charges visible by the electrons, 
the results are close to those obtained by configuration interaction or explicit 
dependence of the wave function on the distance between electrons where tens 
parameters are used. It means that the electron-electron interaction changes 
the electron motion of independent particle mainly between different constant 
interaction potential surfaces affecting the movement along the surfaces to a 
lesser degree. Nevertheless, to improve the accuracy it is needed, in one or 
another way, to take into account the changes in electron movement along the 
surfaces. 

Considering eigenvalues of Q as a functional of parameters determining </), 
it seems possible to avoid the solution of non-linear HF equation when finding 
the best (j) by solving iterative Eq. ([9]). Hopefully, Eq. ([9]) can also replace 
the HF equations for system with more then two electrons. The main effort 
in the proposed technique application is the integration over constant interac- 
tion potential surface. For a two electron system it can be done analytical. 
It doesn't seem difficult to solve this problem for a several electron system. 
However, computation efforts will grow with the increasing number of electrons 
and, eventually, the direct integration in pH)) will be inefficient. We will put 
off discussion the possible ways to solve this problem after solving it for several 
electron systems. 
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TABLE I. Energies of Heliumlike ions (a.u.): HE - HE approximation; Z - 
Eq. ([9]) with effective nuclear charges equal to Zi = Z2 7^ Z, Z\ ^ Z^, ^ Z\ 
CI - configuration interaction. 



ion 






2^5*1 


ion 






2-'Si 


H- 


HE 








HE 


-32.36137° 






Z 


-0.49843 






Z 


-32.38345 


-21.40351 




Z\ = Z2 


-0.50936 






Zi = Z2 


-32.39605 






Zx+Z2 


-0.52623 






Zi + Z2 


-32.40493 


-21.41659 




CI 


-0.52760^ 






CI 


-32.40544'' 


-21.42249= 


He 


HE 


-2.86171'^ 






HE 


-44.73618° 






Z 


-2.87940 


-2.15491 




Z 


-44.75873 


-29.34061 




Z\ = Z2 


-2.89142 






Z\ = Z2 


-44.77422 






Zl^Z2 


-2.90208 


-2.17087 




Zx + Z2 


-44.78001 


-29.35355 




CI 


-2.90325'' 


-2.17521== 




CI 


-44.78061'' 


-29.36118= 




HE 


-7.23633'' 




06+ 


HE 


-59.11159° 






Z 


-7.25642 


-5.09213 




Z 


-59.13392 


-38.52770 




Zi = Z2 


-7.26856 






Z\ = Z2 


-59.15135 






Zl^Z2 


-7.27829 


-5.10636 




Zl^Z2 


-59.15543 


-38.54054 




CI 


-7.27928'' 


-5.11075= 




CI 


-59.15574'' 


-38.55098= 




HE 


-13.61130'^ 






HE 


-75.48702° 






Z 


-13.63244 


-9.27928 




Z 


-75.50906 


-48.96479 




Zi = Z2 


-13.64393 






Zi = Z2 


-75.52834 






Zl^Z2 


-13.65406 


-9.29290 




Zl^Z2 


-75.53066 


-48.97754 




CI 


-13.65485'' 


-9.29739= 




CI 


-75.53083'' 


-48.99223= 




HE 


-21.98607'^ 






HE 


-93.86174° 






z 


-22.00805 


-14.71640 




z 


-93.88415 


-60.65187 






-22.02112 






Zi = Z2 


-93.90520 






Zl^Z2 


-22.02934 


-14.72969 




Zl^Z2 


-93.90577 


-60.66457 




CI 


-22.03020'' 


-14.73463= 




CI 


-93.90592'' 


-60.68527= 



°Ref. [28] 
''Ref. [29] 
=Ref. [12] 
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FIG 1. Efficient charges for l-'^S'o states. 
FIG 2. Efficient charges for 2^51 states. 
FIG 3. Functions x for l^^o and 2^51 states. 

FIG 4. Distribution of the electron-electron distance (solid and dash lines for 
the HF approximation and present calculation, correspondently) and correlation 
holes for 1^5o states. 

FIG 5. Distribution of the electron-electron distance (solid and dash lines for 
the HF approximation and present calculation, correspondently) and correlation 
holes for 2^51 states. 

FIG 6. Function x of the ground state and excited states for non- interacting 
particles. 
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